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Abstract 

We consider the prospects for probing low-mass dark matter with the Super-Kamiokande exper- 
iment. We show that upcoming analyses including fully-contained events with sensitivity to dark 
matter masses from 5 to 10 GeV can test the dark matter interpretation of the DAMA/LIBRA 
signal. We consider prospects of this analysis for two light dark matter candidates: neutralinos 
and WIMPless dark matter. 
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I. INTRODUCTION 



The DAMA/LIBRA experiment has seen, with 8.2o" significance [ij, an annual modula- 
tion [2] in the rate of scattering events, which could be consistent with dark matter-nucleon 
scattering. Much of the region of dark matter parameter space that is favored by DAM A 
is excluded by null results from other direct detection experiments, including CRESST j^, 
CDMS 3, XENONIO j5|], TEXONO (1,0, and CoGeNT 0. On the other hand, astrophysi- 
cal uncertainties j§, [lo| and detector effects 11] may open up regions that can simultaneously 
accommodate the results from DAMA and these other experiments. 



A. A New Window for Light Dark Matter 

The sensitivity of direct detection experiments is largely determined by the kinematics 
of non-relativistic elastic nuclear scattering. If the nucleus recoil energy is measured, the 
momentum transfer can then determined, which in turn determines the dark matter mass 
(for a standard dark matter velocity distribution). This connection allows a direct detection 
experiment to convert an observed event rate into a detection of dark matter with a particular 
mass, or alternatively, convert a lack of a signal into an exclusion of a region of {rriDM, o"si) 
parameter space. 

But the procedure described above implicitly contains the uncertainties which could po- 
tentially lead to consistency between the dark matter interpretation of the DAMA/LIBRA 
signal and the lack of signal at other detection experiments. All direct detection experiments 
have a recoil energy threshold; nuclei recoiling with energies below this threshold cannot be 
detected. Low nuclei recoil energies imply a low dark matter mass. Thus, the sensitivity of 
direct detection experiments tends to suffer at low mass, as seen if Fig. [H DAMA has a 
recoil energy threshold which is lower than many other direct detection experiments, thus 
raising the possibility that its signal is a result of elastic scatter with a light dark matter 
candidate. 

However, there are some direct detection experiments, such as TEXONO [g], 0] and 
CoGeNT [SI], which also have very low energy thresholds, and whose negative results would 
seem to exclude the light dark matter region which would be preferred by the DAMA signal. 
But there are more uncertainties at issue. One example is the measurement of the recoil 
energy. DAMA uses a Nal crystalline scintillator which converts the recoil energy of the 
nuclei into light, which is measured. But, a fraction of the recoil energy is transferred 
to phonons in the lattice and lost. So the measured energy must be scaled by a quenching 
factor to determine the actual nucleus recoil energy. But the DAMA experiment has recently 
noted the existence of the "channeling effect," a property of crystalline scintillators wherein 
a ion moving through certain channels in the lattice will lose no energy to phonons. In 
this case, scaling by the old quenching factor would overestimate the recoil energy, and 
thus the dark matter mass. The DAMA collaboration has adjusted its analysis to account 
for this effect, and their resulting preferred parameter region is not excluded by any other 
experiment [ll|, . 

On the other hand, any uncertainty in the dark matter velocity distribution would alter 
the dark matter mass which should be inferred from any particular measured recoil energy 
distribution. If local streams of dark matter alter the velocity distribution from that expected 
in our local neighborhood, then this could also potentially shift the mass region preferred 
by the DAMA recoil energy signal 
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These effects are detailed in Fig. [H The apphcation of both these effects to the DAMA 
signal is somewhat controversial. It is not at all clear that a dark matter stream with 
sufficient velocity to shift the DAMA preferred region appreciably is reasonable from the 
point of view of astrophysics. And the effect of channeling on a scintillator like DAMA for 
the 1 — 10 GeV energy range is currently being studied and cross-checked by various groups. 

As if this controversy were not enough, three groups 13, 14, 1^ have analyzed the 



spectrum of modulations within the recoil energy bins which DAMA reports. Naturally, 
they all reach different conclusions, ranging from ruling out the dark matter interpretation 
of DAMA to declaring it completely consistent with the spectrum of the DAMA signal. 

It is far from clear what DAMA is actually seeing. What is clear, however, is that if 
DAMA is seeing dark matter, one preferred region of parameter space has dark matter 
mass in the range mx ~ 1 — 10 GeV and spin-independent proton scattering cross section 
o"si ~ 10~^ — 10^^ pb. This is a mass region where several different experimental effects can 
push in different directions, and potentially create a window where dark matter could be 
observed at DAMA while not being ruled out by other experiments. 

Moreover, there are a variety of theoretical models which attempt to explore this region 
of parameter space. Although neutralinos have been proposed as an explanation [l6|, such 
low masses and high cross sections are not typical of weakly-interacting massive particles 
(WIMPs), and alternative candidates have been suggested to explain the DAMA signal [13, 



18|, lia, 120, 121|, 122, 1231, 12J, m 



B. Cross-checking DAMA 

The current state of affairs also makes it abundantly clear that complementary experi- 
ments are likely required to sort out the true nature of this result. Other direct detection 
experiments may play this role. In this work, we note that corroborating evidence may come 
from a very different source, namely, from the indirect detection of dark matter at Super- 
Kamiokande (Super-K) ^26(]. In contrast to direct detection experiments, which rapidly lose 
sensitivity at low masses, Super-K's limits remain strong for low masses. But in contrast to 
other indirect detection experiments, which can only be compared to DAMA after making 
astrophysical assumptions which are highly uncertain (such as the cuspiness of the dark 
matter density profile near the galactic center), Super-K offers a way of testing the DAMA 
result which is largely model independent. Super-K is therefore poised as one of the most 
promising experiments to either corroborate or exclude many dark matter interpretations of 
the DAM A/LIBRA data. 

In Sec. mi we show the relation between the DAMA and Super-K event rates. In Sec. IIIH 
we show that there is significant potential for Super-K to extend its reach to dark matter 
masses from 5 to 20 GeV and provide sensitivity that is competitive with, or possibly much 
better than, direct detection experiments. In Sec. HV] we apply our analysis to two specific 
dark matter candidates that have been proposed to explain DAMA: neutralinos jl6|] and 
WIMPless dark matter 18, 19, 2^- We present our conclusions in Sec. IVl 



II. BOUNDING asi WITH SUPER-K AMIOKANDE 

Super-K can probe dark matter in the Sun or Earth's core annihilating to standard model 
(SM) particles, which subsequently emit neutrinos. Muon neutrinos then interact weakly at 
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or near the detector to produce muons, which are detected at Super-K. The observed rate of 
upward-going muon events places an upper bound on the annihilation rate of dark matter 
in the Sun or the Earth's core. For low-mass dark matter, the dominant contribution to 



neutrino production via dark matter annihilation is from the Sun [27|] , on which we focus. 
The total annihilation rate is 



-Ctanh^ [(aC)^r] 



(1) 



where C is the capture rate, r ~ 4.5 Gyr is the age of the solar system, and a = 
{o'ann.v) / {4:^/2V) , with (Tann. the total dark matter annihilation cross section and V the 



effective volume of WIMPs in the Sun {V = 5.7 x 1030cm3(l GeV/mx)^/^) |27|, l28|, l29|. If 
{o'ann.v) ~ 10~^^ cm^ s~^ (to get the observed dark matter relic density), then for the range 
of parameters considered here, the Sun is in equilibrium {13, [sO, IsH and T 
evaporation is not relevant if mx ^ 4 GeV j2^, [ 
The dark matter capture rate is 27 



30 



^ |C. WIMP 



C 



3n 



Px ^Mb 
a V 

nix 'm 



3(^ 
2 



72/3 



(2) 



The first bracketed factor counts the rate of dark matter-nucleus interactions: a is the dark 
matter- nucleus scattering cross section, px/f^x is the local dark matter number density, m is 
the mass of the nucleus, and Mb is the mass of the capturing object. The velocity dispersion 
of the dark matter is v, and (u^) is the squared escape velocity averaged throughout the Sun. 
The second bracketed expression is the "focusing" factor that accounts for the likelihood 
that a scattering event will cause the dark matter particle to be captured. The parameters /2 
and /s are computable C(l) suppression factors that account for the motion of the Sun and 
the mismatch between X and nucleus masses, respectively, /a ~ 1 for solar capture (27| . 
The capture rate is thus a completely computable function of a/mx- Assuming px = 



300 



km 3 {v^) 



20 [271, and taking /a ~ /g ~ 1, one finds C 



0.3 GeV cm 

10^9 {a/mx) GeV pb"^ s-^ 

The major remaining particle physics uncertainty is the neutrino spectrum that arises 
from dark matter annihilation. Assuming the dark matter annihilates only to SM particles, 
a conservative estimate for neutrino production may be obtained by assuming that the 
annihilation to SM particles is dominated by hb production for < mx < M^, by rf 
production for ruw < mx < mt, and by W, Z production for mx > mt j32[ |. 

Super-K bounds the z/^-flux from dark matter annihilation in the Sun. Since the total 
annihilation rate is equal to the capture rate, this permits Super-K to bound the dark 
matter- nucleon scattering cross section using Eq. ([2]). Fig. [1] shows the published bounds 
from Super-K, limits from other dark matter direct detection experiments and the regions 
of {mx, o"si) parameter space favored by the DAMA signal given astrophysical and detector 
uncertainties. As evident from Fig. [H the published Super-K bounds (solid line) do not test 
the DAMA-favored regions. But we will see that consideration of the full Super-K event 
sample provides significant improvement and extends Super-K's sensitivity to low masses 
and the DAMA-favored regions. 



4 



10 



-2 



CDMS (Si) 



10 



-3 



b 



10-4 



10 



-5 




10-6 



- CRESST 

- TEXONO 

CoGeNT 



WIMPless 



, Super— K 
|v current 



xe non\ \™ ^, (G^) 



5 10 20 
mx (GeV) 




Super-K 
projected 



100 



FIG. 1: Direct detection cross sections for spin-independent X-nucleon scattering as a function 
of dark matter mass mx- The black solid line is the published Super-K exclusion limit 33], and 
the black dashed line is our projection of future Super-K sensitivity. The magenta shaded region 
is DAMA-favored given channeling and no streams [ij], and the medium green shaded region is 
DAMA- favored at 3a given streams but no channeling [l^]. The light yellow shaded region is 
excluded by the direct detection experiments indicated. The dark blue cross-hatched region is the 
prediction for the neutralino models [l^ and the light blue cross-hatched region is the parameter 
space of WIMPless models with connector quark mass my = 400 GeV and 0.3 < Af, < 1.0. Other 
limits come from the Baksan and MACRO experiments 33, 
as Super-K. 
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35| , though they are not as sensitive 



III. PROJECTION OF SUPER-K SENSITIVITY 

As shown in Fig. [H Super-K currently reports dark matter bounds only down to mx 



18 GeV. For heavier dark matter, it is estimated [33[ that more than 90% of the upward- 
going muons will be through-going (i.e., will be created outside the inner detector and will 
pass all the way through it). However, one can study dark matter at lower masses by using 
stopping, partially contained, or fully contained muons, that is, upward-going muons that 
stop within the detector, begin within the detector, or both^. 

Our strategy for projecting dark matter bounds from these event topologies is as follows: 
we begin by conservatively assuming that the measured neutrino spectrum at low energies 
matches the predicted atmospheric background. In any given bin with neutrino events, 
the 2a bound on the number of neutrinos from dark matter annihilation is then 2\/N. This 
bounds the dark matter annihilation rate to neutrinos, which, for a conservative choice of 
the neutrino spectrum, implies a bound on the capture rate, and thus on a /mx- To include 



^ Testing the dark matter interpretation of DAMA has with Super-K through-going muon data has also 
been considered 0. 
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experimental acceptances and efficiencies, we scale our results to the Super-K published 
bound at mx = 18 GeV, assuming these effects do not vary greatly in extrapolating down 
to the 5 — 10 GeV range of interest. 

The annihilation of dark matter particles X typically produces neutrinos with E^/rrix ~ 
I — |. The muons produced by weak interactions of z/^ lie in a cone around the direction to 



the Sun with rms half-angle of approximately 6 = 20°y 10 GeV / [36|. Bounds on dark 
matter with mx = 18 GeV were set using neutrinos with energies ~ 6 — 9 GeV [33| . 
The event sample used consisted of 81 upward through-going muons within a 22° angle of 
the Sun collected from 1679 live days. 

For masses mx ~ 5 — 10 GeV, the z/^ are typically produced with energies between 



2 — 4 GeV. At these energies the detected events are dominantly fully-contained events [37 
so we use this event topology with only events within the required cone around the Sun. 
The number of events is 

,.1 — cos^^ . . 

Nsolar = N , (3) 

where is the total number of fully-contained muon events expected in the 2—4 GeV energy 
range and 9 is the cone o pen ing angle. Super-K expects Ngoiar = 168 such fully contained 
events per 1000 live days [37|. 

We convert this limit on event rate to a limit on the neutrino flux by dividing by the 
effective cross section for the Super-K experiment in the relevant energy range. The effective 
cross-section can be estimated by dividing the estimated rate of events by the predicted 



atmospheric flux, integrated over the relevant range of energies [37[. For fully contained 
events with ii^,^ ~ 2 — 4 GeV, the effective cross section is ~ 2.1 x 10~^ m^. For upward 
through-going events with ^ 8 GeV, the effective cross-section is ~ 1.7 x 10"*^ m^. 

Assuming the neutrino events are detected primarily in either the fully-contained (2 — 
4 GeV) or through-going sample (~ 8 GeV), one can set 2a limits on the time- integrated 
neutrino flux due to dark matter annihilation: 



= , ~ 1.6 X 10^ m-2 

2.1 X 10-8 m2 V 1679 



= , ~ 1.0 X 10^ m-2 , (4) 

1.7 X 10-8 m2 V 1679 ' ^ ^ 

where A^pc = 168 (Adays/1000) and A'tg = 81 (A'days/1679) are the number of fully-contained 
and through-going events within the angle and energy ranges, respectively, scaled to A^ays 
live days. 

The ratio of these flux limits obtained from the fully-contained and through-going samples 
are then equal to the ratio of cr/mx in the 5 — 10 GeV regime to the same quantity at 18 GeV. 
We find 

1.6 X 10^ m-2 / a5_io \ f (^is \ 



1.0 X 109 m-2 \ m5_io 



LSG^vj ' 



where cxs-io is the Super-K bound on the dark matter nucleon cross-section for a dark matter 
particle with mass in the range 5 — 10 GeV, and ais is the bound for a dark matter particle 
with mass 18 GeV. In Fig. [T] this projected Super-K bound is plotted, assuming 3000 live 
days of the SK I-III run. This bound gets better at lower energies and may beat other direct 
detection experiments. 
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IV. PROSPECTS FOR VARIOUS DARK MATTER CANDIDATES 



We now consider specific examples of tlieoretical models that have been proposed to 
explain the DAMA result. We first consider neutralino dark matter. Although neutralinos 
typically have larger masses and lower cross sections than required to explain the DAMA 
signal, special choices of supersymmetry parameters may yield values in the DAMA-favored 
region |16| . 

The region of the (mx, o"si) plane spanned by these models [3] that do not violate known 
constraints is given in Fig. [H We see that if Super-K's limits can be extended to lower mass, 
it could find evidence for models in this class. Note, however, that many of these models 
have p < 0.3 GeV cm~^, and for these models Super-K's bound on the cross section will be 
less sensitive. 

WIMPless dark matter provides an alternative explanation of the DAMA/LIBRA sig- 
nal [l^. These candidates are hidden sector particles that naturally have the correct relic 
density [l8| . In these models, the dark matter particle X couples to SM quarks via exchange 
of a particle Y that is similar to a 4th generation quark. The Lagrangian for this interaction 
is 

C = XfXYLfL + XfXYRfn. (6) 

The Yukawa couplings A/ are model-dependent, and it is assumed that only the coupling to 
3rd generation quarks is significant, while the others are Cabbibo-suppressed.^ One finds 
that the dominant nuclear coupling of WIMPless dark matter is to gluons via a loop of 
6-quarks (t-quark loops are suppressed by m^). The X-nucleus cross section is given by [l^ 



m 



0"si 



N 



2 n2 



E [ZBl + {A-Z)B- 

my — mx 
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(7) 



where Z and A are the atomic number and mass of the target nucleus N , and 5^ 



{2/27)mpfP^/mb, where /|'" ~ 0.8 |38|,l39|. 

In Fig. [l], we plot the parameter space for WIMPless models with my = 400 GeV and 
0.3 < \b < 1.0. These models span a large range in the {mx, ctsi) plane, and overlap much of 
the DAMA-favored region. We see that Super-K's projected sensitivity may be sufficient to 
discover a signal that corroborates DAMA's. But WIMPless models illustrate an important 
caveat to the analysis above; if there are hidden decay channels, then the annihilation rate 
to SM particles is only a fraction of Ttot, and Super-K's sensitivity is reduced accordingly. 



V. SUMMARY 

The DAMA/LIBRA signal has focussed attention on the possibility of light dark matter, 
and alternative methods for corroborating or excluding a dark matter interpretation are 
desired. We have shown that Super-K, through its search for dark matter annihilation to 
neutrinos, has promising prospects for testing DAMA at low mass. 

Using fully contained muon events, we expect that current super-K bounds may be 
extended down to Mom ~ 5 — 10 GeV, and can test light dark models (such as neu- 



tralino models [l6| and WIMPless models ll9|). We have the intriguing prospect that the 



This is a reasonable assumption and is consistent with smaU observed flavor-changing neutral currents. 
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DAMA/LIBRA signal could be sharply tested by an indirect detection experiment in the 
near future. 
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